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DETERMINING A BOUNDARY COEFFICIENT IN A 
DISSIPATIVE WAVE EQUATION: UNIQUENESS AND 
DIRECTIONAL LIPSCHITZ STABILITY 

KAIS AMMARI AND MOURAD CHOULLI 


Abstract. We are concerned with the problem of determining the damping 
boundary coefficient appearing in a dissipative wave equation from a single 
boundary measurement. We prove that the uniqueness holds at the origin 
provided that the initial condition is appropriately chosen. We show that the 
choice of the initial condition leading to uniqueness is related to a fine version 
of unique continuation property for elliptic operators. We also establish a 
Lipschitz directional stability estimate at the origin, which is obtained by a 
linearization process. 
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1. Introduction 

Let SI be a C^-smooth bounded domain of R d with boundary T. We assume 
that r can be partitioned into two disjoint closed subsets with non empty interior 
that are denoted by To and Ti. 

We set, where r > 0 is fixed, 

Q = x (0,t), E o = Pox(0,t), Ei=Pix(0,t), 

and we consider the following initial-boundary value problem (abbreviated to IBVP 
in the sequel) for the wave equation: 

dfu — A u = 0 in Q, 

u = 0 on Eo, 

d„u + b(x)d t u = 0 on Ei, 

u(-, 0) = u°, d t u(-,0) = v°. 

We are interested in the inverse problem of determining the boundary coefficient 
b from the boundary measurement d^Ub^, where Ub is the solution, if it exists, 
of the IBVP (1). At least formally, if each term in the left hand side of the third 
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equation of (1) belongs to L 2 (£i), then the inverse problem under consideration 
is equivalent to determining b from bdtUb^. Therefore this problem is highly non 
linear. 


Before stating our uniqueness and stability results, we need to reformulate the 
IBVP (1) as an abstract Cauchy problem. To this purpose, we set 

v = {w e H x (n); W|r 0 = 0}. 

Here W|r 0 is to be understood in the usual trace sense. When it is equipped with the 
iJ 1 -norm, V is a Hilbert space. We note in addition that the Poincare inequality 
holds true for V and therefore w K > ||Vw|| i 2 ( f2 \ t i defines an equivalent norm on V. 
For s > 0 and 1 < r < oo , we introduce the vector space 

B Str {] R d_1 ) := {w e (1 + \Z\ 2 ) s/2 w e L r (R d - 1 )}, 


where y'(R n 1 ) is the space of temperate distributions on 1 and w is the Fourier 
transform of w. Equipped with its natural norm 

IMIss.rtR^- 1 ) := 11(1 + IC| 2 ) S ^ 2u; llL r (R< 1 - 1 ); 

H Sjr (K d_1 ) is a Banach space (it is noted that B St 2 (M d_1 ) is merely the usual Sobolev 
space i7 s (K" -1 )). By using local charts and a partition of unity, we construct 
B s , r (ri) from H SjT .(R d_1 ) similarly as H s (Ti) is built from H s (R d ~ 1 ). 

The main interest in this spaces is that the multiplication by a function from 
H Sj i(Fi) defines a bounded operator on H s (Ti) (we refer to [3][Theorem 2.1, page 
605] for more details). We set 

B i + / 2 ,i( r i) = {&€ Bi/ 2 ,i(ri); 0 <b}. 


Let b £ 1 (Ti). We define on H = V x L 2 (H) the unbounded operator Ab 


as 


follows 


with domain 


Ah 


0 I 
A 0 


D(Ab) = U v ^ S [H 2 (n)nv] x V; d„u + bv = 0on Tij . 

From [8] [Proposition 3.9.2, page 109], Ab is an m-dissipative operator and there¬ 
fore it generates on T~L a Co-semigroup of contractions e tAb . In that case, we have, 
for any integer k > 1, 


( 2 ) 




u 


v 


0 , 6 f)C k -*([0,T];D(Ai)), 

3=0 

with the convention that A® = I and D(Ab) = %■ 
Moreover, we have the estimate 
,o 


e D(A k b ), 


( 3 ) 


D t^Ab 


V 


< C 


D(A k h ) 


n - =0 cx-iao, t];D(AD) 

Here the constant C doesn’t depend on u° and v°. 

Let A : D(A) C L 2 (fi) -A L 2 (H) be the unbounded operator given by 

A = A, D(A) ={w£ H 2 {9) n V: d u w\ Tl = 0}. 
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The following observation will be useful in the sequel: for any integer k > 1 and 
b £ 1 (Ti), we have 

(4) D(A k ) x {0} C D(A k b ). 


We are now ready to state our main results. 

Theorem 1.1. Let b £ B^ 2 1 (Ti). 

(a) J' = | ^ ^ ^ ; u° £ D(A) and d t uo 7 ^ 0 a.e. on Ei j- 7 ^ 0. 

/ u o \ 

(b) Let I q 1 £ . Then d^Ub^ = d^uo^ implies b = 0. 

Theorem 1.2. We fix b £ B+j 2 non identically equal to zero. We assume 

(5) ( “0 ) G f| D ( A U) 

' ' 0<P<1 

and 

(6) bd t u 0 ± 0. 

Then there exists 0 < po < 1 so that 

n\\pb - 0||s 1/2il (r 1 ) < || d v u p b - 9 1/ mo||l 2( Si ), 0 < p < po¬ 
llere k is a constant independent on p. 

Remark 1. From the proof of Theorem 1.1 (a), we deduce 


jr°o = 


^ ; u° £ and d t uo 7 ^ 0 a.e. on £ 11 7 ^ 0 . 


Therefore, we can replace in Theorem 1.2 (5) and ( 6 ) by 



The authors have already obtained in [1] a log-type stability estimate for the 
inverse problem consisting in determining both the potential and the damping co¬ 
efficient in a dissipative wave equation from boundary measurements. These mea¬ 
surements correspond to all possible choices of the initial condition. The proofs in 
[ 1 ] are essentially based on observability inequalities for exactly controllable systems 
and spectral decompositions. 

The problem of determining a potential in a wave equation from the so-called 
Dirichlet-to-Neumann was map studied by many authors. This problem was initi¬ 
ated by Rakesh and Symes [ 6 ] . We refer the reader who want to learn more on this 
problem to [ 2 ] and references therein. 

The rest of this text is devoted to the proof of our main results. We prove 
Theorem 1.1 in Section 2 and Theorem 1.2 in Section 3. 
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2. Proof of Theorem 1.1 

We first prove a preliminary result. Henceforth, C k denotes the fc-dimensional 
Lebesgue measure. 

Lemma 2.1. Let Asl and u £ H 2 (fi) n C^fl) satisfying 

A u + Xu = 0 in f! and d„u = 0 onT i. 


Then 


C d x ({a; £ ri; u(x) = 0}) = 0. 


Sketch of the proof. Since fl is C^-smooth. Tj can covered by a finite number of 
open subsets U, where U is such there exists a C°°-diffeomorphism if : U —► H, 
B = B( 0,1), satisfying 

i/j(u nn)cH+, ^([/nr)cBo, 

with 

B + = {x = (a/, x d ) £ H; x d > 0}, B 0 = {x = ( x ', x d ) G H; x d = 0}. 

We set v(y) = y G B + . Then Pv = 0 in B + and ddV = 0 on B 0 . Here 

P is a second order operator with C°° coefficients. We extend v to the whole of 
B by setting w = v in B + and w(x',—Xd) = v(x',Xd), (x',Xd) G B + . We have 
w G fl C' 1 (H) because ddV = 0 on Ho and Qw = 0 in B, where Q is a 

second order operator whose coefficients are obtained by taking the even extension 
of the coefficients of P. Checking the details of this construction, we see that Q has 
Lipschitz coefficients. 

Let e > 0 be given. If we denote by TL k the fc-dimensional Hausdorff measure, 
we get by applying [7] [Theorem 2, page 342] that 

H d ~ 2+£ {{z G B- w{z) = 0, \7w{z) = 0}) = 0. 

Therefore, 

n d - 2+e ({y G H 0 ; v(y) = 0, Vv(y) = 0}) = 0, 

In particular 

C. d -\{y G H 0 ; v(y) = 0, Vu(y) = 0}) = 0. 

On the other hand by [4][Lemma 7.7, page 152], V y >v(-, 0) = 0 a.e. in any set where 
n(-,0) is constant. Hence 

C d ~ l {{y G H 0 ; v(y) = 0, d d v{y) = 0}) = 0. 

Bearing in mind that that a Lipschitz map preserves Lebesgue sets of zero measure, 
we get 

C d ~ 1 ({x gPi; u(x) = 0, d v u(x) = 0}) = 0. 

But, d u u = 0 G Ti. Hence 

C d ~ 1 {{x £ Pi; u(x) = 0}) = 0. 


□ 


The following regularity theorem will be useful in the sequel. Since this result 
is not explicitly recorded in the literature, for sake of completeness we sketch its 
proof. 
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Theorem 2.2. Let m > 0 be an integer. For any f £ g £ H m+3 / 2 (To) 

and h £ H m+1 ^ 2 (Ti), the boundary value problem 

{ —Am = / inQ, 

u = g in To, 

d v u = h inT\. 

has a unique u £ H m+2 (Q) satisfying 

(V |l U llTr 2 +™(n) < C (||/||ffm(Q) + || I? || JJm+3/2 (To) + ||ft||ljm+ 1 /2(r 1 )) ■ 

Here the constant C is independent on f, g and h. 


Proof. Since there exists E £ H m+2 (fl ) such that E = g on r 0 and d v E = h in T 
with 

ll- E lli/ 2 + m (fi) <C(\\g\\ H m+ 3 /2(r 0 ) + ||/l||if m +l/2( ri )) ■ 

(e.g. for instance [5][Theorem 8.3, page 39]), we see, replacing u by u — E and / by 
/ — A E, that it is enough to prove the theorem with (g,h) = (0,0). We consider 
then the BVP 

( -Am = / in Q, 

(8) < m = 0 inr 0 , 

I d v u = 0 in F-i. 


Let / £ L 2 {Tl). Because w £ V —> ||Vw||x, 2 (n\d defines an equivalent norm on V, 
by the Lax-Milgram lemma, there exists a unique variational solution u £ H 1 (f2) 
of the BVP (8). That is 


/ Vu-Vvdx= / fvdx 

Jn Jn 


for all v £ V. 


Let (f2o, Hi) be an open covering of an open neighborhood of Q such that C fh 
and r, fl Oi_, = 0, * = 0,1. Let {ipo, tfi) be a partition of unity subordinate to the 
covering (fio, f2i) with ipi £ Cg°(Tli) and ifi = 1 in an neighborhood of Tj, i = 1,2. 

Let Mi = wifi, i = 0,1. Then a straightforward computation shows that Mo and 
Mi are the respective solutions of the variational problems 


/ Vm 0 • T7vdx = / fovdx for all v £ Hq (fl), 
J n Jn 

/ S7u\-S7vdx= j fivdx forall v £ H 1 ^). 
Jn Jn 


Here 


fi = -2V^i- Vu + ipif, i = 0,l. 

Since the regularity theorem [5] [Theorem 5.4, page 165] is valid for both the Dirichlet 
and the Neumann BVP’s, we obtain that m £ H 2 (Tl) and 

IKII/f 2 (fi) < C'll/i||L 2 (n) < C (||M|| ff i(n) + ||/||i 2 (n)) < C"'|l/IU 2 (n), * = 0,1. 

Therefore u = mq + Mi £ H 2 (Tl) and 


IMIff 2 (n) < C'||/|| L 2 ( f2 p 

Next, if f £ 7J 1 (H) then fi £ iJ 1 (H), * = 0,1. We can then repeat the previous 
argument to conclude that u £ 1L 3 (H) and estimate (7) holds with m = 1. We 
complete the proof by using an induction argument in m. □ 
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Proof of Theorem 1.1. (a) Let (A&) be the sequence of eigenvalues, counted accord¬ 
ing to their multiplicity, of the unbounded operator —A. Let (tpk) be a sequence of 
eigenfunctions forming an orthonormal basis of L 2 (f2), each (p k corresponds to A k . 
We note that, according to Theorem 2.2, 

¥>*€ f| H m (Q) = C°°(U). 

m£ N 


We fix k and we take 


Tk 

0 


Then 


G D(A 0 ) ("I D(Ab) and 


uq = cos(i/A ~kt)pk is the solution of the IBVP (1) corresponding to this particular 

/ „,o 

choice of 


We have dtuo = — \/Afc sin(\/A k t)tp k . Hence dtuo ^ 0 a.e. on Id as an immediate 


consequence of Lemma 2.1. In other words, 


Vk 

0 


(b) Let (u°,v°) G J?. If d v Ub | El = S^UoiEj. = 0, then bdtUb^ = 0. Also, by the 
uniqueness of the solution of the IBVP (1), we conclude that Ub = uq. Consequently, 


(9) bd t u 0 | S i = 0. 

By our assumption the set where <9tito|Si vanishes is of zero measure. Therefore, 
(9) implies that 6 = 0 a.e. on Pi. □ 


3. Proof of Theorem 1.2 

We begin by proving an extension lemma. 

Lemma 3.1. (Extension Lemma) Let k. £ two non negative integers. For any g G 
C f ([0, t]; lL fc+1 / 2 (ri)), there exists G G C £ ([0, r]; H k+2 (Ll)) so that, for any t G 
[0,r], 

( AG(t) = 0 in Q, 

< G(t ) =0 on To, 

[ d v G{t)=g{t) onV-L. 

and 

(10) ||GW(t)|U, +2(Q) < C\\g^(t)\\ Hk+1/2{ri) , 0 < j < 1. 

Here the constant C is independent on g. 

Proof. Let h G H k+1 / 2 (Ti). By Theorem 2.2 there exists a unique solution Eh G 
H k+2 (LL) of the BVP 

{ Aw = 0 in Q, 

w = 0 on To, 

d u w = h on Ti. 

Moreover, we have the estimate 

(11) \\Eh\\ H k+2(Q) < C||/l|| ff M-l/2( ri ), 

for some constant C independent on h. 

If g G C^([0,r];IL fe+1 / 2 (ri)), then, using that E : h G H k+1 / 2 ( Ti) -A Eh G 
H k+2 {VL) is linear bounded operator (the fact that E is bounded is a consequence of 
estimate (11)), it is straightforward to check that G(t) = Eg(t) satisfies the required 
properties. □ 
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Next, we consider the following non homogenous IBVP 

( 12 ) 


d 2 u — A u = 0 

in Q, 

u = 0 

on E 0 

d„u = g 

on Ei 

u(-, 0) = u°, dtu(-,0) = v°. 



Proposition 1. We assume that g G C 3 ([0, r]; H l / 2 {Pi)), u° G H 2 (fl) fl V, v° G V 
and the compatihility condition 

(13) d v u° — </(■, 0) = 0 onTi 

holds. Then the IBVP (12) has unique solution u such that 

^ G & = C 1 ([0,t],H 1 (SI) x L 2 (fl))nC([0,T];U 2 (H) x H\n)) 


and 

(14) 


< C 




i7 2 (n)xff 1 (f2) 


+ II5II C 3 ([0,r] ;77 1 / 2 (r i)) ■ 


Moreover, under the additional assumptions 

(15) g G C 6 ([0,T];7J 1 /2 (ri)); ^ ^ 1 ^( 0 ) ) £ D ^' 

v! G C 3 ([0, r]; iJ 1 (J7)) and 

( 16 ) ||^ / ||c' 3 ([0,r];ir 1 (O)) < C ( [ “i ] 

\ \ v J D(A„) 


llffllc £i ([0,T];ff 1 /2(r 1 )) • 


Proof. We denote by G G C 3 ([0, r]; H 2 (VL)) the function given by Lemma 3.1 and 
corresponding to g. We observe that if u is the solution of the IBVP (12) then, 
v = u — G is the solution of following one 


(17) 


Here 


dfv — Av = F in Q , 

v = 0 on So, 

d„v = 0 on Si, 

i>(-,0)= it 1 , d t u(-, 0) = v 1 . 


F = G" 


u 1 = u° - G(0), v 1 =v°- G'(0). 

By the regularity assumptions on u°, v° and g and compatibility condition (13), we 
get that 

F G G 1 ([0, r] ; L 2 (fl)), 'j G D(A„). 

Therefore, the IBVP (17) has a unique solution v so that 

!!, ) eC 1 ([0,r],?{)nC(l0,T];D(A„)). 


This solution is given by 
(18) 


v(t) 

At) 


= e tAo 


M "( 


0 

F(s) 


ds. 
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In light of estimate (10), we have 
(19) 


C 1 ([0,r],W)nC([0,r];_D(A o )) 

,0 


< C 


Since u = v + G, we deduce that 


+ llffllc 3 ([0,r];H 1 /2(r 1 )) I • 


U 

u! 


H 2 (Q)xH 1 (n) 

£ and (19) implies (14). 


Next, we assume that the additional assumptions: 

9 6C 6 ([0,r];H 1 / 2 (r 1 ), ( ) £D(A 4 0 ), 


hold. Then we deduce from (18) that v! £ C 3 ([0, r]; 7J 1 (fl)) and (16) is satisfied. □ 


Proof of Theorem 1.2. We make the following assumption 

)o°)c n d { p ^ 

' 0<p<l 

According to regularity result (2), we have 

u' pb]ri £C e ([0,T];H 1 / 2 (T 1 )), 0 < p < 1. 

and 


( 20 ) 


ll M pb|r 1 ||c 6 ([ 0 ,r];ffi/ 2 (r 1 )) < C 


D(A 7 o) 


, 0 < p < 1. 


We see that v p = u p b — uq solves the IBVP (12) with g = —pbu' pb . By using (16), 
we get 


(21) ll w pllc ,3 ([ 0 > r];.ff 1 / 2 (r 1 )) — ll u pb “ u ollc 3 ([ 0 ,T];fl’ 1 / 2 (r 1 )) < Cp, 0 < p < 1. 

Next, let wo be the solution of the IBVP (12) corresponding to u° = v° = 0 and 
g = —bu'o- Then z = u pb — uq — pwo is the solution of the IBVP(12) corresponding 
to u° = v° = 0 and g = —pb(u' pb — u' 0 ). Hence 

ll^ r llc ,3 ([0,r];ff 1 /2(r 1 )) = P\\K u ' pb ~ M o)llc 3 ([0,T];ib 1 /2(ri))- 
This estimate, in combination with (21), yields 

ll^i'2||c s ([o,T];fl’ 1 /2(r 1 )) < Cp 2 , 0 < p < po- 

Therefore, 

lim — vUpb — U ° = —bd t uo in C 3 ([0, r]; H 1 // 2 (Ti)) 

Pio p 

and then 

hm dvUpb ~ duU ° = -bd t uo in L 2 (Ib). 

Pio p 

By using 2k = ||&d t u 0 || L 2 (El) ^ 0, we get 

Kp < \\d v u p b - ^u 0 ||i 2 ( Sl ), 0 < p < po, 

for some 0 < p 0 < 1 . 

We can rewrite this estimate as follows 

K||p 6 - 0 ||bi /2 i(ri) < \\d v u p b - 9 !/ mo||l 2 (Ei), 0 < p < p 0 . 

This completes the proof. □ 
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